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We study the radiative iieat transfer and tiie van der Waals friction between graphene and an 
amorphous Si02 substrate. We study the surface phonon-polaritons contribution to the low-field 
mobility as a function of temperature and of carrier density. We find that the electric current 
saturate at a high electric field, in agreement with experiment. The saturation current depends 
weakly on the temperature, which we attribute to the "quantum" friction between the graphene 
carriers and the substrate optical phonons. We calculate the frictional drag between two graphene 
sheets caused by van der Waals friction, and find that this drag can induce a high enough voltage 
which can be easily measured experimentally. We find that for nonsuspended graphene the near-field 
radiative heat transfer, and the heat transfer due to direct phononic coupling, are of the same order 
of magnitude at low electric field. The phononic contribution to the heat transfer dominates at high 
field. For large separation between graphene and the substrate the heat transfer is dominated by 
the near-field radiative heat transfer. 

PACS: 47.61.-k, 44.40. -l-a, 68.35.Af 

Introduction. Graphene, the recently isolated single-layer carbon sheet, consist of carbon atoms closely packed in 
a flat two-dimensional crystal lattice. The unique electronic and mechanical properties of graphene^-iS is being actively 
explored both theoretically and experimentally because of its importance for fundamental physics, and for possible 
technological applications^. In particular, a great deal of attention has been devoted to the applications of graphene 
for electronics and sensoringii^. 

Graphene, as for all media, is surrounded by a fluctuating electromagnetic field due to the thermal and quantum 
fiuctuations of the current density. Outside the bodies this fluctuating electromagnetic field exists partly in the form 
of propagating electromagnetic waves and partly in the form of evanescent waves. The theory of the fiuctuating 
electromagnetic field was developed by Ryto\«^"— . A great variety of phenomena such as Casimir-Lifshitz forces^, 
near-field radiative heat transfer^-, non-contact friction^rJ^, and the frictional drag in low-dimensional systems^iii, 
can be described using this theory. 

In the present paper we apply the theories of the near-field radiative heat transfer and the van der Waals friction 
to study the electric transport properties, and the heat generation and dissipation in graphene, mediated by the 
fiuctuating electromagnetic field. The charge carriers in graphene absorbed on, or located close to, a substrate 
experience a friction due to interaction with substrate. This friction affects the graphene carrier mobility and leads 
to the heat generation. The existing microscopic theorieai^ii^ involve several fitting parameters, while our theory is 
macroscopic. In this theory the electromagnetic interaction between graphene and a substrate is described by the 
dielectric functions of the materials which can be accurately determined from experiment. 

The development of more powerful microelectronic devices results in increasing dissipated power (per unit volume) , 
and developing structures with smaller heat production and better cooling strategies, is one of the "grand challenge" of 
modern electronics^^. The heat generated in microelectronic devices leads to elevated device operation temperatures, 
performance reduction, and ultimately to hardware failures. The three classical mechanisms of heat transfer are 
convection, conduction, and radiation. The former is not important in existing circuit architecture having no fiuid or 
gas flowing inside. 

Graphene interact very weakly with most substrates mainly via van der Waals forces. According to theoretical 
calculations^^ the thermal contact conductance due to the direct phononic coupling for the interface between graphene 
and a perfectly smooth (amorphous) Si02 substrate is ctph w 3 x lO^Wm^^K"^, and according to experiment^ (at 
room temperature) the thermal contact conductance ranges from 8x10^ to 1.7x 10^Wm~^K~^ (however, these values 
are probably influenced by the substrate surface roughness). 

At large separation d ^ Xt — ksT/h the radiative heat transfer is determined by the Stefan-Boltzman law, 
according to which a — 4aT^. In this limiting case the heat transfer between two bodies is determined by the 
propagating electromagnetic waves radiated by the bodies, and does not depend on the separation d. At T — 300K 
this law predicts the (very small) thermal contact conductance, a « 6 Wm^^K^^. However, as was first predicted 
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FIG. 1: Two bodies moving relative to each other will experience van der Waals friction due to Doppler shift of the electro- 
magnetic waves emitted by them. 



theoretically^, and recently confirmed experimentall y^°'^^ , at short separation d ^ At, the heat transfer may increase 
by many orders of magnitude due to the evanescent electromagnetic waves; this is often referred to as photon tunneling. 
Particularly strong enhancement occurs if the surfaces of the bodies can support localized surface modes such as surface 
plasmon-polaritons, surface phonon-polaritons, or adsorbate vibrational mode»iSi22,. In two recent experiments22,!^ 
the thermal contact conductance between a silica plate and a silica sphere was measured from 30 nm separation out 
to few micrometers. Heat transfer across the plate-sphere gap causes the cantilever to bent very slightly, and this was 
measured by optical fiber interferometry. 

The theory of the radiative heat transfer developed in Refi^ is only valid for bodies at rest. A more general theory 
of the radiative heat transfer between moving bodies, with arbitrary relative velocities, was developed in KeiM. This 
theory can be applied to calculate the radiative heat transfer between carriers (moving with the drift velocity v) in 
graphene and the substrate. According to theoryi^, for the graphene-Si02 interface the maximum of the contribution 
of the radiative heat transfer to the thermal contact conductivity is Umax ~ {kBT)^n/h, where n is the concentration 



of carriers in graphene. At room temperature and n « 10 cm this gives a 



lO^Wm-^K- 



However this 



estimation does not take into account that the charge carriers in graphene are moving relative to the substrate. 
Relative motion of charge carriers in graphene corresponds to an effective increase in the temperature difference 
between graphene and substrate, what leads to an increase in thermal contact conductance. Thus, for nonsuspended 
graphene the thermal contact conductivity due to near-field radiative heat transfer, and due to direct phononic 
mechanism, can be of the same order. However for suspended graphene, the main contribution results from radiative 
heat transfer since this process exhibit a much weaker distance dependence than for the phononic contribution. 

The radiative heat transfer is closely connected with the van der Waals friction. This friction is due to the Doppler 
effect. Assume that a graphene sheet is separated from the substrate by a sufficiently wide insulator gap, which 
prevents particles from tunneling across it. If the charge carriers inside graphene move with velocity v relative to 
the substrate, a frictional stress will act on them. This frictional stress is related to an asymmetry of the reflection 
amplitude along the direction of motion; see Fig. [TJ If the substrate emits radiation, then in the rest reference frame 
of charge carriers in graphene these waves are Doppler shifted which will result in different reflection amplitudes. The 
same is true for radiation emitted by moving charge carriers in graphene. The exchange of "Doppler-shifted-photons" 
will result in momentum transfer between graphene and substrate, which is the origin of the van der Waals friction. 

Theory. Let us consider graphene and a substrate with flat parallel surfaces at separation d <^ Xt = ch/kBT. 
Assume that the free charge carriers in graphene move with the velocity w <C c (c is the light velocity) relative to 
other medium. According to Refi^rJ^ the frictional stress acting on charge carriers in graphene, and the radiative 
heat flux Sz across the surface of substrate, both mediated by a fluctuating electromagnetic field, are determined by 
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where ni(a;) = [ex^i^'^ I^BTi — 1]^^ (i — g,d), Tgi^^-^ is the temperature of graphene (substrate), Ri is the reflection 
amphtude for surface i for p -polarized electromagnetic waves, and = w ± qxV. The reflection amplitude for 
graphene (substrate) is determined hy^ 

^a{d) = 7~TTT' 
eg(<i) + -L 

where eg(d) is the dielectric function for graphene (substrate). 

In the study below we used the dielectric function of graphene, which was calculated recently within the random- 
phase approximation (RPA)^i2i. The small (and constant) value of the graphene Wigner-Seitz radius rg indicates that 
it is a weakly interacting system for all carries densities, making the RPA an excellent approximation for graphene 
(RPA is asymptotically exact in the <S 1 limit). The dielectric function is an analytical function in the upper 
half-space of the complex w-plane: 

Akpe^ e^g \^( i^ + 2vFkF \ ( oj~2vFkF\ . \ 
eg{uj,q) = l + - —-^====iG\ ]-'^\ ]^^'^(' 

where 

G{x) = x\/l - x2 - ln(.T + v^l - a;2), (5) 

where the Fermi wave vector kp — (TrnY^^, n is the concentration of charge carriers, the Fermi energy ep — "fkp — 
hvpkp, 7 = "hvp ~ 6.5 eVA, and vp is the Fermi velocity. 

The dielectric function of amorphous Si02 can be described using an oscillator model^^ 



where parameters ojqj, and a-j were obtained by fitting the actual e for Si02 to the above equation, and are given 
by too = 2.0014, cTi = 4.4767 x lO^^s^^ wo,i = 8.6732 x lO^s"!, 71 = 3.3026 x lO^^g-i^ = 2.3584 x lO^^s-^, 
a;o,2 = 2.0219 x lO^-^s-^ and 71 = 8.3983 x lO^^s^^ 

The equilibrium or steady state temperature can be obtained from the condition that the heat power generated by 
friction must be equal to the heat transfer across the substrate surface 

F^Td, Tg)v = S,{Td, Tg) + aph{Tg - T^), (7) 

where the second term in Eq. ([7]) takes into account the heat transfer through direct phononic coupling; Uph is the 
thermal contact conductance due to phononic coupling. 

Results. At small velocities [v <^ vp) the friction force depends linearly on v: = Tv, where the friction 
coefficient F is determined by 
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FIG. 2: The role of the interaction between phonon polaritons in Si02 and free carriers in graphene for graphene field-effect 
transistor transport. The separation between graphene and Si02 is d = 3.5A. (a) Dependence of the low-field mobility on the 
temperature. Concentration of charge carriers n = lO^^cm"^^. (b) Dependence of the low- field mobility on the concentration 
of charge carriers n for Td^Tg = 300 K. (c) Current density-electric field dependence at T = K, n = lO^^cm"^^ (d) The 
same as in (c) but for different temperatures. 



The low-field mobility can be described using Matthiessen's rule 

f^^^ = f^int + t^ext, (9) 

where fiint is the intrinsic mobihty due to scattering of the graphene free carriers against the acoustic and optical 
phonons in graphene, and iiext is the extrinsic mobility due to scattering against the optical phonons in the substrate. 
According to the theory of the van der Waals friction fiext = ne/T. Scattering of the graphene carriers by the acoustic 
phonons of graphene places an intrinsic limits on the low-field room temperature mobility given hy fi — 200000cm^/Vs 
at the carriers density lO^^cm^^ However, for graphene on Si02 the mobility of the graphene carriers is reduced 
to fj, — 40000cm^/Vs. This reduction must result partly from interface imperfections, and partly from scattering of 
the graphene carriers from the optical phonons of the Si02 substrate. 

Fig. [2^ shows the calculated dependence of the low-field mobility on the temperature, due to scattering from 
the optical phonons in Si02. The calculations were performed using the theory of the van der Waals friction. The 
mobihty is characterized by a strong temperature dependence. At room temperature the calculated mobility 
3 X 10^ cm^/Vs, which is approximately one order of magnitude larger than the mobility calculated in Ref.^^. This 
indicates that for graphene adsorbed on Si02 the main contribution to the low-field mobility comes not from scattering 
from the optical phonons in Si02, but from scattering from point defects on the interface. 

Fig. shows the dependence of the mobility at Td — 300 K on the charge density n. The mobility monotonically 
decreases when the concentration of charge carriers increases. Figs. [2]: and[2ji show the dependence of the current 
density on the electric field at n = lO^^cm"^, and for different temperatures. In obtaining these curves we have used 
that J = nev and neE — , where J and E are current density and electric field, respectively. Note that the current 
density saturate at E ^ 0.5 — 2.0V/ /im, which is in agreement with experiment^. The saturation velocity can be 
extracted from the I — E characteristics using Jsat — nev sat, where 1.6 mA//im is the saturated current density, and 
with the charge density concentration n = lO^^cm^^: Vsat ~ lO^m/s. Fig. ^ was calculated at = K. At zero 
temperature the van der Waals friction is due to quantum fluctuations of charge density, and is determined by the 
second term in Eq. ^-22^ 

F^ATd = = 0) = --3 / dq, \ dq., / dojq^e-^'^'^ ..^^o^..-M. (10) 
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FIG. 3: A schematic diagram of a drag experiment. A current Ji is passed through layer 1, and a voltmeter is attached to 
layer 2 to measure the induced electric field E2 due to the inter-layer interaction. 

The existence of "quantum" friction is still debated in the literature^"—. The theory^rJ^i^ predicts that a solid moving 
relative to another experiences a force due to quantum fluctuations, that is opposite to the direction of motion. The 
van der Waals friction can be studied in non-contact experiments, and in frictional drag experiment a 1 ^ ^ . In both 
these experiments the solids are separated by a potential barrier thick enough to prevent electrons or other particle 
with a finite rest mass from tunneling across it, but allowing the interaction via the long-range electromagnetic field, 
which is always present in the gap between bodies. In non-contact friction experiments the damping of cantilever 
vibrations is typically measured, while in frictional drag experiments a current density is induced in one medium. The 
friction between the moving charge carriers and nearby medium gives rise to a change oi I — E characteristics, which 
can be measured. 

The friction force acting on the charge carriers in graphene for high electric field is determined by the interaction 
with the optical phonons of the graphene, and with the optical phonons of the substrate. The frequency of optical 
phonons in graphene are a factor 4 larger than for the optical phonon in Si02. Thus, one can expect that for graphene 
on Si02 the high- field I — E characteristics will be determined by excitations of optical phonons in Si02. According 
to the theory of the van der Waals frictio n '^"i^^ , the "quantum" friction, which exist even at zero temperature, is 
determined by the creation of excitations in each of the interacting medium, the frequencies of which are connected 
by vqx = LJi-\-uj2- The relevant excitations in graphene are the electron- hole pairs which frequencies begining from zero, 
while for Si02 the frequency of surface phonon polaritons ujph ~ 60meV. The characteristic wave vector of graphene is 
determined by Fermi wave vector kp. Thus the friction force is strongly enhanced when v > Vsat = '^ph/kp ^ lO^m/s, 
in the accordance with numerical calculations. Thus the measurements of the current density-electric field relation of 
graphene adsorbed on Si02 give the possibility to detect "quantum" friction, what has fundamental significance for 
physics. 

An alternative method of studying of the van der Waals friction consists in driving an electric current in one metallic 
layer and studying of the effect of the frictional drag on the electrons in a second (parallel) metallic layer (Fig. |3]). 
Such experiments were first suggested by Pogrebinskii'^^ and Price—, and were performed for 2D-quantum wells^i^. 
In these experiments, two quantum wells are separated by a dielectric layer thick enough to prevent electrons from 
tunneling across it, but allowing inter-layer interaction between them. A current density Ji = niev is driven through 
layer 1 (where ni is the carrier concentration per unit area in the first layer), see Fig. [3] Due to the inter-layer 
interactions a frictional stress a = Tv will act on the electrons in layer 2 from layer 1, which will induce a current 
in layer 2. If layer 2 is an open circuit, an electric field E2 will develop in the layer whose infiuence cancels the 
frictional stress a between the layers. Experiments^ show that, at least for small separations, the frictional drag 
can be explained by the interaction between the electrons in the different layers via the fiuctuating Coulomb field. 
However, for large inter-layer separation the frictional drag is dominated by phonon exchange'^^. 

Similar to 2D-quantum wells in semiconductors, frictional drag experiments can be performed (even more easy) 
between graphene sheets. Such experiments can be performed in vacuum where the contribution from phonon exchange 
can be excluded. To exclude noise (due to presence of dielectric) the frictional drag experiments between quantum wells 
were performed at very low temperature (T w 3 K)2i. For suspended graphene sheets there are no such problem and 
frictional drag experiment can be performed at room temperature. In addition, 2D-quantum wells in semiconductors 
have very low Fermi energy ep ~ 4.8 x 10~'^eV— . Thus electrons in these quantum wells are degenerate only for very 
low temperatures T < Tp = 57 K. For graphene the Fermi energy ep — O.lleV at 71 = lO^^cm"^, and the electron 
gas remains degenerate for T < 1335 K. 

At small velocities the electric field induced by frictional drag depends linear on the velocity, E — {T /ne)v = n^^v, 
where fi is the low-field mobility. For huj <^ ep and q ^ kp the reflection amplitude for graphene is given by the same 
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FIG. 4: Frictional drag between two graphene sheets at the carrier concentration n = lO^'^cm"^. (a) Dependence of friction 
coefficient per unit charge, ^T^ = F/ne, on the separation between graphene sheets d. (b) Dependence of electric field induced 
in graphene on drift velocity of charge carriers in other graphene sheet at the layer separation d = 1 nm. (c) The same as in 
(b) but at d = lOnm. 



expression as for a 2D-quantuni welliS; 



Substitution of Eq. (HIl) into Eq. (g]) gives 
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Fig. shows the dependence of the friction coefficient (per unit charge) on the separation d between the sheets. 
For example, E = 5x 10~'^v for T — 300 K and d = 10 nm. For a graphene sheet of length 1 fj,m, and with v — lOOm/s 
this electric field will induce the voltage y = 10 nV. Figs. [SId and[5l; shows the induced electric field- velocity relation 
for high velocity, with d = Inm (b) and d — lOnm (c). 

As discussed above, for graphene on Si02 the excess heat generated by the current is transferred to the substrate 
through the near-field radiative heat transfer, and via the direct phononic coupling (for which the thermal contact 
conductance a « lO^Wm-^K"^). At smaU temperature difference (AT = Tg -Td<^ Td), from Eq. Q we get 



where F^o = F^{Td,Tg = Td), S^o = Sz{Td,Tg Td), 
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The thermal contact conductance is given by 

_ S^Td^Tg) + gp/iAT _ {aph + S'^Q)Fr,QV ~ SmF'j,qV 



AT 



Fxqv — Szo 



(14) 




FIG. 5: Radiative heat transfer between graphene and Si02. (a) The dependence of the ratio between the total heat generated 
by current and the radiative heat flux, on electric field, n — lO^^cm"^, d = 3.5A, = 1.0 x 10*Wm~^K~^ (b) Dependence 
of graphene temperature Tg on the electric field. The substrate temperature Td = 300 K. Parameters are the same as in (a), 
(c) Dependence of the thermal contact conductance on the separation d for low electric field, aph ~ 0. 



From Eq. (fT4|) it follows that the thermal contact conductance can be strongly enhanced when Fxov w Szo- Fig. [5^ 
shows the ratio of the total heat generated by the current to the radiative heat flux. For low field this ratio ~ 2 so 
that in this case both radiative heat transfer and phononic heat transfer give contributions of the same order. The 
phononic mechanism of heat transfer dominates for high electric field. For high electric field the equilibrium graphene 
temperature due to self-heating depends linear on the electric field (see Fig. [Sja) in contrast to the low-field case 
where this dependence is quadratic. This can be explains by the fact that the saturation velocity depends weakly 
on electric field; therefore heat, generated by friction, will linearly depend on the electric field. Fig. [S}; shows the 
dependence of the thermal contact conductance on the separation d for low electric field. At d 5 nm and T = 300 
K the thermal contact conductance, due to the near-field radiative heat transfer, is ^ 10'*W-m^^-K^^, which is ^ 3 
orders of magnitude larger than that of the black-body radiation. In comparison, the near-field radiative contact 
conductance in Si02-Si02 for the plate-plate configuration, when extracted from experimental data^*' for the plate- 
sphere configuration, is ^ 2230W-m^^-K^^ at a ^-^30 nm gap. For this system the thermal contact conductance 
depends on separation as l/cP; thus a ~ 10^W-m^^-K~^ at d '-^ 5 nm what is one order of magnitude larger than 
for the graphene-Si02 system in the same configuration. However, the sphere has a characteristic roughness of ~ 40 
nm, and the experiments^!^ were restricted to separation wider than 30 nm (at smaller separation the imperfections 
affect the measured heat transfer). Thus the extreme near-field-separation, with d less than approximately 10 nm, 
may not be accessible using a plate-sphere geometry. On the other hand suspended graphene sheet has a roughness 
~1 nm^I, and measurements of the thermal conductance can be performed from separation larger than 1 nm. At 
such separation one would expect the emergency of nonlocal and nonlinear effects. This range is of great interest for 
the design of nanoscale devices, as modern nanostructures are considerably smaller than 10 nm and are separated in 
some cases by only a few Angstroms. Another advantage of using graphene for studying the radiative heat transfer 
result from the fact that under steady-state condition the heat fiow is equal to the heat generated in the graphene by 
the current density: Sz — EJ. This quantity can be accurately obtained from I — V characteristics. The temperature 
of graphene can be measured accurately using Raman scattering spectroscopy'^^. 

Conclusion. We have used theories of the near-field radiative heat transfer and the van der Waals friction to 
study transport, heat generation and dissipation in graphene due to the interaction with phonon-polaritons in the 
(amorphous) Si02 substrate. In contrast with existing theories, based on the semiclassical Boltzmann transport 
equation, our approach is macroscopic. In the latter approach all microscopic properties are included in dielectric 
functions of material. Explicit formulas were obtained for the low-field mobility, the high-field transport and the 



near-field radiative heat transfer between graphene and the substrate. The low-field mobility exhibit a strong tem- 
perature dependence, which (according to the theory of the van der Waals friction) is associated with the thermal 
fluctuation inside the media. High-field transport exhibit a weak temperature dependence, which can be considered 
as manifestation of quantum fluctuations. Thus the study of transport properties in graphene gives the possibility to 
detect "quantum" friction, the existence of which is still debated in literature. We have calculated the frictional drag 
between graphene sheets mediated by the van der Waals friction, and found that it can induce large enough voltage 
to be easily measured experimentally. This effect can be used in electronics and for sensoring. We have shown that 
for the low-field heat transfer between graphene and the substrate, both radiative heat transfer and phononic heat 
transfer give contributions of the same order. High-field heat transfer is determined by the phononic mechanism. We 
have pointed out that graphene can be used to study near-field radiative heat transfer in the plate-plate configuration, 
and for shorter separations than it is possible now in the plate-sphere configuration. 
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